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ABSTRACT 

In a quadruply imaged lens system the angular distribution of images around the 
lens center is completely described by three relative angles. We show empirically that 
in the three dimensional space of these angles, spanning 180° x 180° x 90°, quads from 
simple two-fold symmetric lenses of arbitrary radial density profile and arbitrary radially 
dependent ellipticity or external shear define a nearly invariant two dimensional surface. 
We give a fitting formula for the surface using SIS+elliptical lensing potential. Various 
circularly symmetric mass distributions with shear up to 7 ~ 0.4 deviate from it by 
typically, rms ~ 0.1°, while elliptical mass distributions with ellipticity of up to e ~ 0.4 
deviate from it by rms ~ 1.5°. The existence of a near invariant surface gives a new 
insight into the lensing theory and provides a framework for studying quads. It also 
allows one to gain information about the lens mass distribution from the image positions 
alone, without any recourse to mass modeling. As an illustration, we show that about 
3/4 of observed galaxy-lens quads do not belong to this surface within observational 
error, and so require additional external shear or substructure to be modeled adequately. 



1. Introduction 

Recovering projected mass distribution of galaxies and clusters given the images of lensed 
background sources is an important problem, and much effort has been devoted to lens mass 
modeling over the last couple of decades. In this paper we show that useful information about the 
lensing object can be obtained without any recourse to mass modeling. We work exclusively with 
quadruply imaged lens systems (the fifth central image is usually not detected and is not part of the 
analysis), and more specifically, with the angular distribution of the four point-like images around 
the lens center, that was introduced in Williams et al. ( 2008[ ). We do not consider image fluxes. 



A typical quad image configuration is shown in Figure [TJ The images are labeled by their 
arrival time at the observer, 1 through 4. In most cases this ordering can be determined from 



the morphology of the quad, without measured time delays (Saha Sz Williams 2003). Here we 
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are interested only in the mass distribution of the lens, not its total mass (or, equivalently, its 
normalization), or orientation. In this case the image configuration of any quad is described uniquely 
by 6 parameters, which we chose to be of the polar variety, measured with respect to the lens center: 
three relative angles, and three distance ratios of images. The three relative angles between images 
are marked on the plot, #12, #34 and #23. Angle #12 is between the two minima of the arrival time 
surface, while #34 is the angle between the saddle points. We define #12 such that it encloses image 
3, and 6*34 such that it encloses image 2. Of the three angles, #23 is special because the separation 
between images 2 and 3 gets arbitrarily small for sources approaching the diamond caustic shown in 
the lower right panels of Figure [TJ when the source crosses the caustic these two images disappear 
and the quad becomes a double. Note that any linearly independent combination of the above three 
angles can be used, but we chose 612, #34 and 6*23 because they have a simple physical meaning. 



Working with only two angles, #23 and a certain linear combination of #12 and #34, Williams et 



al. (2008) showed that a wide range of simple, twofold symmetric lens models generate apparently 
indistinguishable patterns in the two dimensional plane of these angles. Twofold symmetric means 
that the mass distribution, and hence the potential, is symmetric about two orthogonal axes, 
and 'simple' excludes lenses with 'wavy' isodensity shapes. The simple, twofold symmetric class of 
lenses includes all popular parametric lens models, such as Singular Isothermal Ellipsoids (SIE), and 
Singular Isothermal Elliptical Potential (SIEP), as well lenses of any density profile and ellipticity. 



The present paper is an extension of Williams et al. ( 2008 ) , but here we work with the full 



set of three angles, #12, #34 and 623- We show that quads from all simple lens mass distributions 
with twofold symmetry lie on nearly the same two dimensional surface in the three dimensional 
space of relative angles. We call this the Fundamental Surface of Quads (FSQ). The quads from 
observed galaxy lenses, on the other hand, show a different behavior. As we show in Section [5j 
galaxy quads form a 'cloud' surrounding the FSQ, with typical separations from the FSQ of few to 
several degrees. 

One can draw some interesting parallels between the FSQ we introduce here and the well 
studied Fundamental Plane of Ellipticals. Both lie in the three dimensional space whose axes are 
parameters describing the structural properties of the respective objects. In the case of quad lenses, 
these are the relative image angles, while in the case of ellipticals they are the effective radius, the 
surface brightness at the effective radius, and the central velocity dispersion. A wide class of objects 
belong to the Surface and the Plane with small scatter. In other words, the objects do not fill the 
full three dimensional space, implying that there is a tight relation between the three parameters. 
The existence of the Fundamental Plane is basically the consequence of the virial theorem, while 
the reason for the near invariance of the Fundamental Surface of Quads for a wide class of twofold 
symmetric lenses (but not necessarily for the observed quads) is yet to be identified. 
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2. The SIS+elliptical lensing model 



We start by studying a simple, analytically tractable, two dimensional projected gravitational 
potential. It belongs to the generic family of separable potentials, <j>(r,0) = r-f(0), where r and 
are polar coordinates in the lens plane. Properties of such potentials are discussed in Kassiola & 



Kovner 


(1993 


), 


Kochanek 


(1991 


) and Dalal|(1998) 



hence, 



/(0) = &[l + 7 cos(20)] 
cj) = rb[l + 7cos(26>)] 



(1) 
(2) 



which is sometimes called SIS+elliptical; we will call it SISell for short. The normalization factor 
b is the Einstein radius, and 7 is the shear parameter. 



The Poisson equation, A(p = 2k, yields the projected dimensionless mass density profile, 

K = — [1-37 COS (20)1. 
2r 



(3) 



Note that 7 cannot be greater than 1/3 since otherwise k will have an unphysical negative value. 
The lens equation, 

r s =r-V(j) (4) 

where f s and r are source and image positions respectively, can be rewritten as two independent 
equations 

r s cos(0 - 9 S ) = r - b[l +7cos(26»)] (5) 



r, sin 



-2&7sin(20) 



-a sin(# — 6 S ) = sin(20), a 



267' 



(6) 



Setting the magnification M = 1/ det(^4) , where A is the Jacobian matrix of the lens equation, 
to infinity, i.e. det(vl) = 0, one gets 



det(A) = - [(r - b{l + 7 cos(26»)}) + 46 7 cos (20)1 = 0, 
r 

which yields the condition for the caustic, 

r = b-3b-/cos{26). 



(7) 



(8) 



Now using equation ^ in the lens equation, eq. [6j allows one to express the caustic coordinates 
(r sc , SC ) in the source plane as a function of parameter 0, 



6 7v / 2[5 + 3cos(40)] 
- tan -1 (tan 3 9) . 



(9) 
(10) 
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Our aim is to calculate the three relative angles, #12, #34 and #23 for all the quads within 
the diamond caustic. Using eq. [6j which is independent of image distance r, we get the angular 
positions of the four images, in radians. 



j -1 / 1 L , -a 3 cos 3 ((9 s ) + a(a 2 -4)cos(# s ) + 8acos(0 s ) 1 5 \ 

1 = COS 2 V AB ^acos{8 s ) + - \ , (11) 



, -1 / 1 -a 3 cos 3 (# s )+a(a 2 -4)cos(# s )+8acos(# s ) 1 5\ 

9 2 = cos I — - y C — -acos(0,)-— I , (12) 

a L -a 3 cos 3 (0 s ) + a (a 2 - 4) cos(^) + 8a cos(0 s ) 1 . s\ 

#3 = cos -acos(0 fl ) -- 1 , (13) 



1 / 1 L, -a 3 cos 3 (6U + a(a 2 -4)cos(0 s ) + 8acos(6U 1 ,„ . S\ , 
#4 = -cos" 1 I -2V C + — — K -^--a C os{9 s ) + -\ , (14) 

where, 



A = f/ 108a 4 sin 2 (20 s ) + 2 (a 2 - 4) 3 + 12\/3Wa 4 sin 2 (2# s ) (27a 4 sin 2 (2fl s ) + (a 2 - 4) 3 ) , (15) 



c - ! w + ? u «^ ( -')-4 (17) 

The above equations are rather cumbersome and so preclude simple analytical expressions for 
the relative angles. Instead, we numerically generate random source positions within the caustic, 
calculate 9i's using the above equations, and then compute relative angles. 

The symmetry of the potential implies that it is sufficient to consider source positions only 
within one of the quadrants of the elliptical potential. Therefore, without any loss of generality we 
align the x-axis with the major axis of the ellipse and consider 9 S only from to ir/2. For a given 
9 S , r s can vary in the range [0,r sc ), therefore with eqs. [6] and [9] we see that a runs in the range of 
0, 0.5-^/2 [5 + 3 cos(4#)] J , which is independent of 6 and 7, where the parameter 9 is determined 



by the value of 9 S using eq. 10 



3. The Fundamental Surface of Quads 



We use the expressions for 9^s given in the previous Section to parametrically plot the relative 
quad angles in the three dimensional space of #12, #34 and #23- The resulting distribution is a two 
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dimensional surface, shown in Figure [2^ a). The fact that it is a surface means that in a quad 
resulting from a SISell lens, two relative image angles completely determine the third. 

The surface is simple with distinct properties. It has a slightly curved triangular shape with 
its apex at (#12, #34, #23) = (180°, 180°, 90°); quads at the apex have a "cross" configuration. The 
two edges connecting the apex to the base at 6*23 = 0° correspond to #12 = 180° and #34 = 180°, 
because in a twofold symmetric lens the latter two angles do not exceed 180°. Note that the base of 
the triangular surface that represents small values of #23 and source positions close to the caustic, 
shows some unevenness, or jaggedness. It is unclear whether this is intrinsic to Q{ equations, or if 
it is due to the numerical noise arising from the implementation of these equations. In all what 
follows we ignore these small features; in particular, our fit to the surface, discussed below, smooths 
over this unevenness. 

This surface is universal for all SISell lenses, however the meaning of universality requires 
some clarification. For a given source position (9 S , r s ), a relative angle 9ij depends on b and 7, and 
therefore two SISell models with different b and 7 give rise to two different points on the surface. 
However, the surface itself does not depend on b and 7. This is the result of the elimination of b 
and 7 dependence when considering all source positions within the caustic as discussed in Section 
[2] Therefore quads from SISell models of all shears and normalizations lie on the same invariant 
surface. 

We would like to have an explicit functional form for the surface, as #23 = f cn (^i2 ; #34)) but 



since the equations for individual angles, eq. 11 14 contain inverse cosines and are complicated 



there is no simple expression. Instead, we calculate thousands of sets of relative angles, (#12, #34, 
#23)) from the expressions for the 0j's (eq. l~lfl~4 ) and fit these with a surface represented by a 



polynomial function in #12 and #34. The fitting was done using Matlab's Least Absolute Errors 
(LAE) method. As compared to the Least Squares method, LAE is less stable and could generate 
more than one function. We chose LAE anyway because it is resistant to outliers, which in our case 
correspond to quads with small values of #23, and are responsible for the jaggedness of the surface. 

We determine the optimal order of the polynomial by considering the deviations, or errors, of 
the SISell quad points from the fit surface, quantified by the root mean square error, RMSE. A 
number of trials and tests revealed that a fourth order polynomial] has the lowest value of RMSE, 
~ 0.00032 radians, or ~ 0.018°, for approximately 160 thousand quads. We note that using different 
sets of quads to obtain the fit equation resulted in slightly different values for the coefficients of the 
polynomial, and for some sets of quads the RMSE varied by up to a factor of two. 

To test the robustness of our fit to changes in the fitting method, we also computed the best 
fit surface using the Least Square method^ and using custom versus in-build polynomials within 



1 Matlab allows up to 5th degree polynomial fit. 

2 Numerical Recipes' Singular Value Decomposition routine did not provide a good fit when single precision was 
used, while in double precision svdcmp failed to invert the matrix at all. 
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Matlab. Again, the resulting fit surface changed somewhat, but did not deviate substantially from 
the fit we present below. We conclude, therefore, that although the fit equation is not reproducible 
exactly, it is completely adequate for our purposes: 

(9 23 = -5.792 + 1.783 6 12 + 0.1648 6 2 l2 - 0.04591 6\ 2 - 0.0001486 9\ 2 + 1.784 34 

-0.7275 34 0i2 + 0.0549 34 0j 2 + 0.01487 34 0? 2 + 0.1643 0f 4 + 0.05493 0f 4 9 12 (18) 
-0.03429 0f 4 - 0.04579 0| 4 + 0.01487 0| 4 i2 - 0.0001593 0| 4 

Figure [2^b) shows the fit as the gray semi-transparent surface. The red points are SISell quads 
corresponding to a random distribution of source positions within the diamond caustic, on the 
source plane. As depicted in the Figure, a random distribution of source positions does not imply a 
random distribution of quads on the Fundamental Surface; quad density increases with increasing 
023- Two other orientations of the Fundamental Surface are shown in Figure [3j 

Because the RMSE of the SISell quad distribution about the fit plane is < 0.02° the differences 
between the two will be invisible in the full three dimensional angles space. Instead, we calculate 
the difference in 2 3 of the SISell quads and the fitted surface keeping the other two angles fixed; we 
call this difference A023 = 023 — &23 jit, where 9 2 3ju is obtained by plugging 0i2 and 03 4 of a quad 
in to eq. 18 Figure [4] plots A023 vs. 2 3. The straight horizontal line represents the surface fit, 



eq. 18, while the points are the SISell quads. The wiggles in the distribution of points represents the 



wiggles in the SISell surface, compared to the fourth order polynomial fit. However, for all practical 



purposes, the differences are negligible, and eq. 18 can be taken to be a good representation of the 
SISell potential. 

We call this surface the Fundamental Surface of Quads because, as we show in the next section, 
not just SISell, but most twofold symmetric models do not differ from it by more than a few degrees. 
This near invariance probably stems from the shape of the caustic of this class of potentials. The 
twofold symmetry of the lensing mass distribution implies the twofold symmetry of the diamond 
caustic. More specifically the diagonals of the caustic intersect at 90° and all four quadrants of the 
caustic are identical. 



4. Other two-fold symmetric potentials 

In this Section we explore a wider range of simple twofold symmetric mass models, motivated 
by the commonly used parametric models. We calculate typical values for a total of 12 models, but 
show plots (see below) for only eight of these. 

Two of the first four models have isothermal (SIS) radial density profiles, and the other two, 
de Vaucouleurs (deV) profile. Isothermal means that if mass ellipticity were zero, the projected 
density profile would scale as 1/r. The de Vaucouleurs profile has projected density given by, 



£ = S e exp I —7.673 



(r/r e ) 



1/4 



(19) 



-7- 



where r e is the half- mass radius, and S e is the projected mass density at r e . To generate a 
diamond caustic the density profiles must be accompanied by either ellipticity, e, or external shear, 
7. Ellipticity, e of the mass isodensity contours is related to the axis ratio of the isodensity contours, 
b/a = (1 — e)/(l + e). The properties of the first four models are: SIS with e = and 7 = 0.3; deV 
with e = and 7 = 0.4; SIE with e = 0.3 and 7 = 0; deV with e = 0.4 and 7 = 0. 

The additional eight models explore a range of power-law surface mass density profiles, oc 
R~ a , with q = 0.4, 0.7, 1.3, 1.6, where isothermal slope is a = 1. The range of slopes we have chosen 
is considerably wider than what real lensing galaxies seem to have. Based on a well defined sample 



of 15 elliptical galaxy lenses, Koopmans et al. (2006) find that the typical total (dark matter and 



baryons) space density slope is 2.01, or about 1 is projection, with dispersion of 0.12. We chose 
a range of slopes 5 times wider than that because we would like to demonstrate the robustness of 
the FSQ to changes in the lens model parameters. Each of these four density profile slope models 
were given e = 0.25, or 7 = 0.25. These values are somewhat on the high side of the typical values 
for ellipticity and shear encountered in modeling observed quads. 

Each mass model generates a two dimensional surface in the three dimensional space of relative 
angles. All surfaces coincide at the apex, where (9\2, #34, #23) = (180°, 180°, 90°), since all twofold 
symmetric models can produce a perfect cross configuration (though the distance ratios of images 
will differ). At other locations the surfaces deviate from each other somewhat, and from the 
one defined by SISell; largest deviations are seen at small #23 values, i.e. towards the base of 
the triangular surface. However, even at the widest separation, the surfaces differ by only a few 
degrees. The left panel of Figure [5] shows that these differences are hard to discern in the full three 
dimensional angles space, even if the surface is split up into four pieces for easier visualization. To 
make the deviations visible, in the right panel of Figure [5] we fold the surface along the vertical 
mid-line, and show only a narrow range of angles, a few degrees in each case. In this zoom, the 
deviations are seen to be a few degrees. 

Another way to show deviations is to use A623 introduced in the previous Section; see Figures [6] 
and [7} As before the horizontal line at A623 = represents the Fundamental Surface of Quads, 
given by the fit eq. 18, and the points are quads from lens mass models. In Figure [6] the isothermal 



models and the de Vaucouleurs models are shown in the left and right panels, respectively. The 
top row shows the models with external shear, 7, while the bottom rows represent elliptical mass 
distributions. In Figure [7] four mass models with power-law density profiles are shown. 

To quantify the effect of shear or ellipticity on any given mass model, we quote the average 
distance ratio < d^/di >, where d\ and d& are the distances of the first and fourth arriving images 
from the lens center, and the average is over quads randomly populating the inside of the diamond 
caustic. A given lens mass distribution produces quads with a range of image distance ratios, but 
in general image 1 tends to be farthest from the lens center, while image 4 tends to be closest. For 
the 12 models, < d^jd\ > is between 0.5 and 0.9. For the sample of 40 observed quads (Section [5]) 
< d^/di >~ 0.68. i.e. typically smaller than in our models. As will be discussed in the next 
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Section many of the observed quads are the result of mass distributions that are more involved 
than two-fold symmetric lenses; most require an external shear in addition to and elliptical lens, 
while some require substructure. The presence of these would tend to reduce the d±/d\ ratio. 

Table [T] summarizes the results of the 12 models. In general, larger ellipticities or larger 
7 result in larger deviations from the FSQ. Because the maximum deviations from FSQ differ 
between models from a fraction of a degree to a few degrees, the range on the vertical axes range 
are different in the top and bottom panels of Figures [6] and [7j As opposed to the quad distributions 
generated by the elliptical mass models, the ones from models with external shear are much closer 
to the FSQ, and appear more similar to that of the SISell model (unless the surface mass density is 
very shallow, as in the bottom right panel of Figure[7|). When viewed in 3D, the surfaces containing 
quads from the elliptical lens models sag below the FSQ, but even for the bottom right panel of 
Figure [6] the deviations in #23 are < 4°, for e = 0.4, and de Vaucoulers profile. Here the images 
are formed where the projected density oc R~ 2 , or p(r) oc r -3 in three dimensions. Given 

observed lenses, this is a rather extreme combination of ellipticity and T,(R). Ellipticity of e = 0.4 
corresponds to the axis ratio b/a = 0.43, or an E5.7 if it were an optical elliptical galaxy. Steep 
density profiles also appear to result in larger deviations from the FSQ; see Table [TJ Real galaxy 
lenses rarely, if at all, have such steep profiles at the location of quad images. For e = 0.25, and 
oc #-0-4...-l.6 (t p right and bottom left panels of Figure [7]), the deviations from FSQ are a 
degree at most. 

We note that for very large ellipticities or shears (not considered here) the two opposite cusps of 
the diamond caustic protrude outside of the oval caustic producing so-called naked cusps, which do 
not produce quads. The corresponding surfaces of relative angles look similar to the ones without 
the naked cusps, except that the portions at the bottom corners of the surface are devoid of quads. 



5. Observed Quads 

In this Section we illustrate one of the practical uses of the Fundamental Surface of Quads 
(FSQ). 

Galaxy lens systems can be approximately divided into three categories, depending on whether 
the lens mass model is (a) an elliptical mass distribution or a circularly symmetric mass distribution 
with an external shear, (b) an elliptical mass distribution plus some external shear, or (c) a more 
complicated mass distribution, possibly with additional lens galaxies or substructure mass clumps. 
A survey of the literature indicates that only a handful of systems belong to (a) . A model- free way 
to come to that conclusion is to look at the quads in the 3D angles space. 

We have assembled a sample of 40 galaxy-lens quads. The sample was collected from all 
available sources, and is therefore heterogeneous. Where possible, the astrometry, including the 
positional uncertainties on the images and the lensing galaxy was taken from the CASTLeS web- 



site (Kochanek et al. 2008); otherwise from individual papers. In the latter case, systems listed in 
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Table [2] have a footnote with a reference; a lens system with no reference means that its data were 
obtained entirely from CASTLeS. The image arrival time was determined from the morphology 



of the lens (Saha & Williams 2003), and the relative angles were calculated. These are listed in 



Table [2j in columns labeled #12, #34 and 623- 

In Figure [8] we show two orientations of the 3D angles space with the FSQ and the 40 quads. 
For clarity, we do not show errors in this plot. The main conclusion is that most observed quads 
lie more than a few degrees away from the FSQ; 12 are within ±2°, so most cannot be modeled 
adequately with an elliptical lens, or a circularly symmetric lenses with external shear. 

Next, we incorporate errors into the analysis. Even though the astrometric measurement errors 
of images and galaxy lens center are largely independent of each other, a shift in the lens center 
translates into correlated relative angle errors. To account for this we calculate the errors as follows. 
We assume that the x, y positional errors of each of the four images and the lens center are normally 
distributed, with a x and <r y taken from the literature. We then draw thousands of independent 
image and lens center positions, and for each generated lens system calculate relative angles #12, 
#34 and #23- The thousands of generated quads per lens system then give us the error distribution 
for each of the three relative angles. We calculate the mean and the rms of these distributions, and 
list them in columns labeled 0ij^ errors in Table [2j Note that the average of these distributions need 
not be the same as 6ij, however, the differences tend to be small, generally < 0.1°. 

We quantify the deviation of the quads from the FSQ as in earlier Sections, by calculating 
A(9 23 . The error in A$23) listed as A$23, errors is calculated similarly to what was described above, 
using thousands of quads generated based on astrometric uncertainty. 

Figure [9] shows the distribution of the 40 quads in the #23 vs. A#23 plane. Within errorbars, 
only 10 systems are consistent with FSQ. Several of these have published parametric modeling, and 
are, in fact, well represented by two-fold symmetric mass distributions. For example, B2045+265 
is successfully modeled by Fassnacht et al. ( 1999 ) using SISell potential, eq. [2] SDSS J002240 is 



modeled by Allam et al. (2007) with Singular Isothermal Ellipsoid (SIE) using gravlens software 



of Keeton (2001). The same software was used by Grillo et al. (2010) to fit the positions (not the 
flux ratios) of SDSS J1538 with three types of twofold symmetric models: de Vaucoulers, SIE and 
a power law density profile. 

On the other hand, some of the lenses which lie away from the FSQ are known to require 
external shear in addition to elliptical lens. PG 1115 where the lensing galaxy is a member of a 
galaxy group is inconsistent with FSQ; its |A6>23| ~ 4°. RXJ 0911+0551 has a cluster next to it 
( Burud et al.fl998 ), so the lens model requires an external shear in addition to an elliptical galaxy 
lens; its |A0 23 | ~ 6°. For LSD Q0047-2808, iKoopmans & Treul (12003b state that SIE+shear does 



not fit the image positions well (but sufficient for the determination of the Einstein ring radius); 
it has |A#23| ~ 14°. Lenses with known secondary galaxies also lie far from the FSQ. HE 0230- 



2130 has a secondary lensing galaxy (Wisotzki et al. 1999) very close to the images, and so its 
|A#23| ~ 30°. B1608 has a complicated galaxy merger as a lens, and its |A#23| ~ 5°, and so it is 
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inconsistent with the FSQ. 

A few caveats are in order. If a quad does not lie within a couple of degrees of the Fundamental 
Surface of Quads (i.e. the range defined by the various density profile and ellipticity models, such 
as the ones in Figure [5j [6] and [7]), it cannot be modeled by a twofold symmetric lens. However, 
the opposite need not be true. If a quad lies on the FSQ does that immediately imply that it can 
be modeled by a twofold symmetric lens, regardless of its image distance ratios? This question 
will be addressed in a future study. We also note that even if a quad does belong to the locus of 
twofold symmetric lens in the full 6D space of image position properties, it does not mean that 
other types of lens models cannot fit it. Reconstruction of the lens mass reconstruction from single 
quads is a highly underconstrained problem, so the solution is not unique, and many mass models 



can reproduce the image positions exactly (Saha & Williams 2004) 



6. Conclusions and Future Work 

In this paper we present a model-free way of making inferences about the lensing mass dis- 
tribution given its quad image positions. The latter are represented by three relative angles that 
describe the distribution of images around the lens center. We show that in the three dimensional 
space of these angles, quads generated by SIS+elliptical mass distribution belong to an invariant 
two dimensional surface, regardless of the shear parameter 7, and normalization b. Furthermore, 
quads from a wider class of lenses with twofold symmetry outline almost the same surface, making 
the surface a near invariant descriptor of twofold symmetric mass distributions. Because of that 
property we call the two dimensional surface the Fundamental Surface of Quads (FSQ). 

The existence of FSQ allows one to characterize galaxies and clusters based on the quads they 
generate. To aid in that, we provide a fitting formula for the FSQ based on the SIS+elliptical 
lensing potential. If a quad does not lie within a couple of degrees of the FSQ (i.e. the range 
defined by the various density profile and ellipticity models, such as the ones in Figure [5] [6j and 
[7}, the mass distribution is not twofold symmetric. This method of determining if a lens can be 
fit with a twofold symmetric lens is superior to answering this question using parametric modeling. 
The latter fits quads with a finite set of models, while our method addresses all twofold symmetric 
models irrespective of the specific parametric form. 

However, the main importance of the FSQ is not in ascertaining if the lens mass distribution is 
twofold symmetric or not, but in the following aspects, which we will investigate in the forthcoming 
papers. First, the near invariance of FSQ provides a new framework for studying quads, and strong 
lensing theory in general. We remind the reader that it is still not understood why a wide class 
of twofold symmetric lenses form such a tight, near invariant distribution in the space of relative 



angles. Second, as already shown in Williams et al. (2008), the relative angles present a promising 



way of investigating realistic mass distributions, and specifically, differentiating substructured lenses 
from smooth non- two fold symmetric ones. Finally, the full set of quad image properties lives in the 
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six dimensional space that includes image distance ratios. An investigation of this space is yet to 
be undertaken. 
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Fig. 1. — A representative quad from a twofold symmetric lens. Top left: An elliptical lens mass 
distribution, with the n = 1 contour shown as a thick blue line. Images are the magenta filled 
circles. Top right: Arrival time contours and images. Bottom left: Images, labeled by arrival time. 
The relative angles, #12, #23 and #34 are marked. Bottom right: The diamond caustic, with the 
location of the source represented by a solid blue square. 



Fig. 2. — The three dimensional space of three image angles, 6u, #23 and #34 for the SIS+elliptical, 
or SISell mass distribution, (a) The surface outlined by quads whose relative image angles were 
calculated using eq. 11 14 (b) The fit surface to the Fundamental Surface of Quads (FSQ), eq. 18 



is shown as the gray surface, while the points are quads randomly distributed within the diamond 
caustic in the source plane. 
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Two additional orientations of the Fundamental Surface of Quads. 
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Fig. 4. — Deviations of the SISell quads from the best fit 4th degree polynomial presented in equa- 
tion 18 The vertical axis shows the deviations of the quads' 023 from the prediction of equation 18 
The quads used to compute the best fit were obtained using analytical equations for angles 6i, 82, 
83 and 84 presented in Section [2j while the quads in this Figure were generated using a ray tracing 
code. The difference between the two is small. 
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Fig. 5. — Simple twofold symmetric lens mass distributions define a nearly invariant Fundamental 
Surface of Quads. Quads from four mass models are shown: SIS with e = and 7 = 0.3 (red); deV 
with e = and 7 = 0.4 (blue); SIE with e = 0.3 and 7 = (black); deV with e = 0.4 and 7 = 
(brown). On the left we show the 3D space of relative angles sliced into four segments divided 
by #23 = 35° , 55° , 75° . The fact that the points of different lens potentials are hard to tell apart 
demonstrates the near invariance of the FSQ. To make the deviations visible, on the right we fold 
and project a small angle range of the surface; see Section [4] for details. 
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Fig. 6. — The deviations of the quads of four mass distributions from the predictions of the 4th 
degree polynomial fit, equation [TFJ The four mass distributions are the same as the ones shown in 
Figure [5] Top left: Circularly symmetric SIS surface mass density with external shear 7 = 0.3; Top 
right: Circularly symmetric de Vaucoulers with external shear 7 = 0.4; Bottom left: Elliptical SIE 
with ellipticity 0.3; Bottom left: Elliptical de Vaucoulers with ellipticity 0.4. The average value of 
the distance ratio of the fourth to first arriving image, < d^/di >, is shown in each panel. This 
aids in visualizing the meaning of 7 and e value. Note that the vertical axes have different ranges 
in the top and bottom panels. 
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Fig. 7. — Same as Figure [6j but for a different set of lens mass models. Top left: Circularly 
symmetric surface mass density, X(-R) oc R with external shear 7 = 0.25; Top right: Elliptical 
E(R) oc R~ 0A with ellipticity 0.25; Bottom left: Elliptical E(R) oc J? -1 - 6 with ellipticity 0.25; 
Bottom left: Circularly symmetric oc R~ 0A with external shear 7 = 0.25. 
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Fig. 8. — The Fundamental Surface of Quads fit equation (shaded region) and the forty observed 
galaxy quads (red dots). Two orientations are shown; in the second one the deviation of the 
observed quads from the Plane are clearly visible. 
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Fig. 9. — Deviations, A#23, of the observed quads from the FSQ. The horizontal axis is the observed 
023- The empty circles represent the observed relative angles and their deviations from the FSQ. 
The error bars are calculated as explained in Section [5} The horizontal and vertical axes values are 
given in Table [TJ 



